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Abstract 

In this paper we show how to compute in a consistent way nucleation rates in 
field theory at finite temperatures, with metastable vacuum. Using the semiclassical 
approach in field theory at finite temperature we show that the prefactor term 
can be calculated explicitly (in the thin-wall approximation) and that the same 
provides exponential finite temperature quantum corrections to nucleation rates, 
when fluctuations around the bubble field configuration are considered. 
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1 Introduction 



The study of first-order phase transitions in cosmology has received increased interest 
recently due to their possible relevance to the physics of the early Universe. Indeed, 
phase transitions have been invoked in the context of the generation of large scale structure 
of the Universe, the flatness and horizon problems' 1 ' and, more recently, the eletroweak 
phase transition, which could be responsible for the generation of the cosmological baryon 
asymmetry^ . 

In a first-order phase transition the Universe is believed to be in a metastable phase 
(the false-vacuum) that decays to a stable phase (the true vacuum) as the Universe ex- 
pands and cools. The decay process of the metastable phase to the stable one happens by 
the nucleation in the system of droplets (or bubbles)^ of the stable phase whose dyna- 
mics will determine the completion of the phase transition' 4 ^. Of fundamental relevance 
for the understanding of the development of the phase transition is the determination of 
the bubble nucleation rate per unit volume. 

In this paper we deal with the problem of determination, using the semiclassical ap- 
proach, of the nucleation rate in field theory at finite temperature, which includes its 
definition from a droplet model field theoretical point of view' 3 ! and the evaluation of the 
prefactor term appearing in the nucleation rate expression. We show that the prefactor 
term provides an exponential finite temperature quantum correction to the nucleation rate 
which has not been properly discussed in most of the studies involving this quantity M'^. 
Our method of evaluation is self-consistent and we show that the prefactor can be eva- 
luated directly by computing the eigenvalues of the determinant terms, in the case of a 
convenient bubble field configuration, or by developing an appropriate field theoretical 
expansion for the determinant ratio. 
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The usual expression for the nucleation rate per unit volume used in the literature is 
given by! 5 ! 




the false vacuum (or metastable) field configuration, S 3 ((f),T) is the three-dimensional 
Euclidean action, at finite temperature, of the non-trivial field configuration 0, 



where V e ff(<p,T) is the effective potential at finite temperature. 

What is usually done to compute the determinant ratio in (1) is the use of dimensional 
analysis to approximate it by a preexponential factor of order T 4 , for example, T 4 (the 
critical temperature) or m 4 (0j,T) (m 2 (0/,T) = d Vef /Jf'' T ^ \<j>=<j> f )- This way of defining 
the nucleation rate presupposes that the determinant ratio in (1) will not contribute to 
the term in the exponential and that all the quantum corrections are included, from the 
beginning, in the form of S^((p,T) in the exponential. 

We shall show that our way of defining the nucleation rate T results in a expression 
quite different of that given by (1) and that the usual approximation for the overall 
preexponential factor is not a simple function of the temperature, especially near T c , 
as observed recently by Csernai and Kapusta^, that have computed the preexponential 
factor in the case of a QCD phase transition. 

The paper is organized in the following way: in Sec. II we review the droplet model 
in field theory at finite temperature for a generic model of a scalar field 0, which could 
be coupled with other fields, gauge or others scalar fields, for example. If the effective 
potential at finite T, of the scalar field 0, admits a first-order phase transition, then 
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the nucleation rate can be defined as being the imaginary part of the free energy, in the 
approximation of a dilute gas, of droplets of the stable phase inside the metastable phase' 3 ' . 
This imaginary free energy arises due to the existence of a negative eigenvalue, associated 
with the instability of the critical bubble. In Sec. Ill we deal with the evaluation of 
the prefactor appearing in the nucleation rate and show that the same is of fundamental 
importance for the proper determination of the finite temperature corrections to the 
exponential factor term in the nucleation rate. In Sec. IV we study the model of a scalar 
field theory coupled with fermion fields, computing the aditional determinant factor due 
to the integration over the fermion fields. In Sec. V we apply the method by choosing 
a simple potential which exhibits metastability, at T ^ 0. Conclusions are presented in 
Sec. VI. 

2 Droplet Model in Field Theory 

Let us consider a system described by a scalar field <p and by a set of fields x (bosonic or 
fermionic fields), that can be coupled to <fi. The partition function Z of this system, at 
T ^ 0, is 

Z = J D(f)Dxe~ SEucl{M , (3) 

where Se U ci((P, x) is the Euclidean action of the system and the functional integrations are 
carried over field configurations subject to periodic (anti-periodic) conditions: <p(x, 0) = 
4>(x, f3) for bosons and if)(x, 0) = —ip(x, (3) for fermions. (3 is the inverse of the temperature, 

If one integrates out the x fields in (3) and makes the expansion <f>(x, r) — > (f) c (x, r) + 
rj(x, r), where </> c (x, r) is a field configuration that extremizes the Euclidean action S Euclid) 
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(after integrating out the x fields) for the scalar field <J) and r](x, r) are small perturbations 
around <f> c (x, t), then the effective action r e //(0 c ) for the field configuration C , can be 
written as 

T e// (0 C ) = In | D V e- §E ^ , 
or using an expansion in terms of <p c and its derivatives, 

r e/ /(0c) = jfV I d 3 * [-y e// (0 c ,r) + i(a M c ) 2 ^(0 c ) + ...]. (4) 

In (4), for constant field configurations, <j> c (x, r) = C , V e ff((f) c ,T) defines the effective 
potential, at finite T. In principle, from V e ff(<p,T) we can determine the order of the 
phase transition. If above a certain critical temperature T c the system is in a symmetric 
phase 4>f and below T c V e ff((f),T) develops an energetically favorable phase <p t , then the 
phase (pf becomes a metastable phase and we have a first-order phase transition in which 
the metastable phase 0/ decays to the new favorable phase (f> t . T c is defined by the 
condition V eff ((f) f ,T c ) = V eff ((f) t ,T c ). 

The process of the decay of the metastable vacuum can be seen as happening by the 
nucleation of bubbles (or droplets) of the new phase 4>t i* 1 the system, that in a successful 
phase transition will grow and coalesce completing the phase transition. 

Let us write the partition function (3), after integrating out the x fields as 

Z = ^D0exp{-5 BttC ,(0)} , (5) 

where Se U ci{.<P) is an effective Euclidean action for the scalar field <f>. 

If the system defined by (5) admits an effective potential like that one described above, 
then one can imagine the transition from the metastable phase <pf to the stable phase <f> t 
as being preceeded by the formation of small nuclei of the phase <f> t inside the phase (j)f. 



Note that <fif and <fi t are functions of the temperature, <fif( t ) = 0/( t )(T), given by the 



minima of V eff (<f>, T), H=4>m 



0. Following Langer^, in a dilute gas of droplets 



of the phase (fit , one can infer the thermodynamics of the system from the knowledge of 
the partition function of a single bubble. Therefore, in a dilute gas approximation, the 
principal contribution for the partition function (4) for the system can be written as 

Z ~ Z((f> f ) + Z(<k) 

and the partition function Z of this gas of bubbles can be approximated to' 3 ' 



z ~ z{4> f ) 



z i<t>f) ex p 



(6) 



where in the expressions above Z(<f>f) is the partition function for the metastable phase 
and Z(<pb) is the partition function for the bubble in the system, or in a field theory 
language, Z(<p b ) is the partition function of the system in the presence of a bubble field 
configuration (fib ' 7 '- 

For S E ud(4>) = Jo dr J (i 3 x[|(5 M 0) 2 + V(4>)], the bubble field configuration is a non- 
trivial solution of the Euler-Lagrange equation: 



Or 2 



+ V 2 )(f>-V'(<fi) = 



(7) 



and therefore, in the semiclassical approach, the field configuration <fi b is the one that 
extremizes the effective action 







(8) 



For static field configurations, from Eq. (7), <fib(r) is a static bounce solution of the 
differential equation 
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with boundary conditions, lim r ^ +00 = </>y and ^| r =o — 0. 

One can then describe 0&(r) by the following spherically symmetric field configuration 

(j>t, 0<r<R-AR 
Mr) = l Kaiu R-AR<r<R + AR (10) 
(p f , r>R + AR 

that is, 4>b(r) describes a bubble of radius R, of the nucleating phase (fi t embedded in the 
metastable phase <fif with the bubble wall described by a field configuration (f) wa ii, that 
separates the two phases (for example, one could imagine 4> wa ii as a kink-like field configu- 
ration). AR is the bubble wall thickness. The solution (10) is a very good approximation 
for fa when R » AR, or in the thin- wall approximation. The simplest assumption 
for the validity of the thin-wall approximation is to consider that the bubble nucleation 
happens at temperatures close to T c , that is, with relatively small supercooling®. In the 
rest of this work we will concentrate on this assumption, where we can approximate the 
bubble field configuration quite well by <fo(r) given by (10). 

Let us now compute the partition function (6) using a semiclassical approach. Within 
the semiclassical approach one expands the Lagrangean field (fi in (6) as (fi — > r)(x, t)+<f>b(x) 
for Z((fib) and (fi — > ((x, t) + (fif for Z(<fif). rj(x, t) and ((x, t) are small perturbations around 
the classical field configurations and <fif respectively. Up to 1-loop order one keeps 

the quadratic terms in the fields rj(x,t) and ((x,t) in these expansions. In this way one 
can write the following expressions for Z(<fi b ) and Z(<fif), respectively, 

Z((fi b ) 1 - lo % order e -s E M ^DrjexpL^dr J d 3 x^r) [-U Eucl + V" (</> b )] rj\ (11) 
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and 



zm ^e-^) /D C exp{-/ o ^r / d*x\t [-D Eucl + V"(<f> f )] C | , (12) 



where S Eud {4>) = fg ' dr J d?x[\{d,4>) 2 + V(<j>)], V"{<t>) = and U Eu(d = + V 2 . 

Performing the functional gaussian integrals in (11) and (12) one gets the following 
expression for the ratio between the partition functions, z(4> b f \ ' a PP earm g in (6): 



Z{^>b) 1 — loop order 



det(-D Eud + V"{<j) b ))p 



(13) 



_det(-a Eucl + V"(<p f )) l3 _ 
where [det(M) i g] - ^ = § Drjexp {- dr / ^x\q[M\n} and AS = S Eucl (<f) b ) - S Eucl ((j> f ) is 
the difference between the Euclidean actions for the field configurations (fi b and 0/ . 

The free energy of the system, T = — \nZ, up to 1-loop approximation, from (6) 
and (13), can be written as 



T = —T 



det(-n Eucl + V"((t) b ))i 



(U) 

det(-n Eucl + v»(<p f ))p\ ■ 1 ; 

It is well known that the determinant for the bubble field configuration in (14) has 
a negative eigenvalue, that signals the presence of a metastable state, and that it has 
also three zero eigenvalues related with the translational invariance of the bubble in the 
three dimensional space. Because of the negative eigenvalue, the free energy T is then 
imaginary. However the imaginary of T can be related exactly with the nucleation rate 
of bubbles of the phase <f> t (the stable vacuum) inside the metastable phase <f>f, as shown 
by Langer' 3 '. At finite temperatures, Affleck, ref. [8], showed that the nucleation rate T 
is given by 



where \cu-\ is the frequency of the unstable mode. 



3 Evaluation of the Determinants 

Let us now compute the ratio of the determinants appearing in the nucleation rate, Eq. 
(14), and show that the same provides a finite temperature correction to the exponential 
term AS . 

One remembers that AS is given byQ 

r A E 

AS = f3 J d 3 x [C Eud {<f) b ) - £ Eucl {<P f )} = — . (16) 

If one uses the bubble field configuration </>&(r) as given by (10), then in the thin- wall 
approximation one obtains the following expression for AE in (16): 

AE = ^-AV + 4nR 2 a , (17) 

where AV = V((j)f) — V(<pt) is the potential difference between the false and true vacua 
(the metastable and stable vacua, respectively) and <7q is the surface tension of the bubble 
wall (with no corrections due fluctuations around the bubble wall field configuration 4> wa ii) 



(To 



r+AR _ 

/ dr C E ud{(t>wall) - £ E ucl((pf) ■ (18) 

J-AR L J 



Remind that in the above equations and in <fib( r ), given by (10), </>/( t ) = 0/( t )(T), are the 

minima of the finite temperature effective potential, that is, <f>f(t) are given by 

1 We are supposing that £,Euci(4>) denotes an "effective Euclidean action" for the field <j>, where possible 
other fields coupled with <fi has been integrated out. In this way V(<fi) includes not only the classical 
potential for cf> but also can include corrections (T ^ 0) coming from the integration of that fields. 
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— ^ — l^/(t)-°- I 19 ) 

From (17) one can see that AE can be associated with the activation energy of one 
bubble of radius R. What we are going to show is that the determinant ratio in (14) will 
give a finite temperature correction to this bubble activation energy, where this correction 
comes exactly from fluctuations around the bubble field configuration c/) b . 

The computation of the ratio of determinants in (14) can be done by using two ap- 
proaches. The first one involves obtaining directly the eigenvalues of the determinants in 
Eq. (14). The second one consists in developing a consistent expansion for the determi- 
nant ratio. 

3.1 Evaluation of the prefactor in terms of the eigenvalues 

The evaluation of the determinants in (14) can be done by computing directly, if possible, 
the eigenvalues of the differential equations 

\-U Eud + V"(<j> f )} <p f (i) = e 2 f (i) Vf (i) (20) 

and 



[-U Eud + V" (fa)] <p b (i) = e? b (i)<p b (i) ■ (21) 
In momentum space one writes, e 2 = u 2 + E 2 , where uj n = , n — 0, ±1, ±2, for 
bosons (for fermion fields uj n = <y2n+ ^ 1T ). From (20) and (21) one can write the determinant 
ratio in (14) as 

' det(-a i ; Md + y // (0 / )) / , 
_det(-0 Eucl + V"((f) b ))p 




det(-a gMd + y // (0/)) /3 
det(-n Eucl + V"((j) b ))p 
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Using the identity: 



exp { - In 



n 



+00 



11 K 



E 2 Ai 



n 



+00 



n, H 



Em 



(22) 



n 1 



sinh(7rz) 



n=l 



(23) 



rr ) Ti z 

and taking into account that we have in (22) a negative and three zero eigenvalues, one 
obtains for (22) the expression: 



det(-a^ + K // (^/)) / 3 
det(-n Eucl + V"((f) b ))p 



^-Isin 



AS 



2tt 



x exp<j]T 
"0 



f^«+ln(l-e-W)) 



^E b (j) + In (l - e -^W) 



(24) 



The factor 



AS 



2tt 



in the right hand side of (24) comes from the contribution of the zero 
eigenvalues ^ . The prime in J2j is a reminder that we have excluded the negative and the 
three zero eigenvalues from the sum. E 2 _ denotes the negative eigenvalue (|u;_| = \E_\, in 
(15) ). 

Substituting (24) in (14) and using the bubble field configuration given by (10), 

one obtains the following expression for the nucleation rate T, given by (15): 



T = „4T 4 exp 



AF(T) 
T 



(25) 



where we have denoted by A the adimensional factor a 



AE 
2vrT 



with a given by 



a = — 



\E-\ 
2T 



(26) 
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and AE is given by (17). AF(T) in (25) is given by 



AF(T) = -^-AU(T)+A7iR 2 a(T) , 



(27) 



where 



d 3 k 



In 



d 3 k 

1 _ e -f3^fP+m?{<t>t) 



AU{T) = V(<P f ) - V{4>t) +Tj In [l - e-^+mH^) 



(28) 



and 



a(T) 



^o + ^|E ' ^E wall {j) + In (l- e 
Eff^)+ln(l-e-W))j| . 

i L J ) 



■/3E wal i(j) 



(29) 



In (27) we have used again the thin-wall approximation. In (28) we have substituted the 
discrete sums by integrals over momenta (for the constant field configurations <pf and (j) t , 
we have the continuum eigenvalues, Ej = k 2 + m 2 {(f)f) and E 2 = k 2 + m 2 {(f> t ) , respectively, 
with m 2 (0 f ) = ^^U= 0/ and m 2 ^) = ^^1^ )• The terms, like / d 3 k^Jk 2 + m 2 (0), 
that should appear in (28) are ultraviolet divergent and can be subtracted from the theory 
by the introduction of the usuals counterterms of renormalization that render the theory 
finite^. 

In (29) E wa u(j) are the eigenvalues related with the bubble wall field configuration 
(f) wa u- The problem then is reduced to the computation of these eigenvalues for a field 
configuration describing the bubble wall. 

It is easy to see that AU(T) as given by (28) is exactly the finite temperature effective 
potential difference between the false and true vacua, as expected^ 11 '. The second term in 
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the right hand side in (29) clearly represents the finite temperature contribution for the 
surface tension o"o, coming from the 1-loop finite temperature quantum corrections due 
to fluctuations around the bubble wall field configuration (j) wa u. 

3.2 A field theoretical expansion for the determinants 

The second approach for the computation of the determinants in (14) consists in develop- 
ing a simple field theoretical expansion for it. Let us write for the ratio of the determinants 
the following expression 



det(-a gMc/ + r / (0/)) /3 

det'(-n Eud + y"(0 6 ))^ 



= exp |-Tr In [-D Eud + V"{4> f )] p - 
- \Tr' ln[-D BucJ + V"(^} , 



(30) 



where we have used in (30) the identity In det M = Tr In M and the prime in both sides 
denote that the negative and the zero modes have been omitted. 
Formally one can write (30) as 



' det(-a i ; Mc/ + r / ((/'/)) / 3 

det'(-a Eucl + V"{<f> h ))p_ 
where 



_ i 

2 



= exp 



-Tr In 



1 + Gp((f)f) [V"(0 6 ) " V"(<f>f)\ 



(31) 



^ -n M + m >(+ f ) (32) 
is just the free propagator, at finite temperature, for the scalar field (f), with mass squared 

given by m 2 (0/) = V"(<f> f ). 

If we expand the natural logarithm in (31) in powers of Gp{(j>f) \Y"{<j>b) — V'i&f)], we 

get formally 
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Tr \n{l + G^ f )[V"( ( p b )-V"( ( f )f )]} = 



+ ... , (33) 



where the dashed lines correspond to the background-like field [V"((pb) — V"{4>f)\ an d the 
internal lines denote the propagator Gp(<f>f). The expression (33) can be written as 



+00 / i \m+l „ 

Tr In {1 + G^ f ) [V"(<f> b ) - V" (</>,)]} = £ { -^— / d 3 x [Vfa) - V"(<P f )] m x 

m=l m J 

^ 1 , , 

X n ^J (2,r)3[ a; 2 + p + m2(0/) i m - ^ 
The sum in m in (34) can be formally performed and one obtains 



Tr ln{l + G /3 (0 / )[^"(0 6 )-y"(0 / )]} = J d'x £j _ln 



1 , - 

+ k 2 + m 2 ((pf) 

' (35) 



Substituting (35) in (31) and again taking into account that we have eliminated the 
negative and zero modes from the determinant in (14), which provides us with that factor 
multipling the exponential in (24), one obtains an expression for the nucleation rate T 
like (25), but now with AF(T) given by (using (10) as a solution for <p b and the thin- wall 
approximation) 



AF(T) 



AV eff (T)+±7rR 2 a(T), 



(36) 
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where 



AV eff (T) = V(4 f ) - vfa) - 



+00 



20 ^ 

iJ t J n=— oc 



d 3 k 

(2tt) 3 



In 



1 + 



m 2 ((j) t ) -m 2 ((P f ) 
u 2 + k 2 + m 2 ((j)f) 



and 



(37) 



a(T) = a + 



1 



1 +00 



In 



1 + 



V"( 



'wall J 



u 2 + k 2 + m 2 ((j)f) 



(38) 



AttR 2 J 20„±n„J (2tt) 3 
Expression (37) can easily be identified with AU(T), given by (28), if one performs 
the sum in n in (37) by using the identity 



+00 



<4 + E?(k) 
" 2 n + E}(k) 



= pE t (k) + 2\n(l-e- f>Et ®) - 
- 0E f (k)-2\n(l-e- f3E f {%) ) 



(39) 



and as it was done in (28), eliminating the ultraviolet divergent terms, we obtain 



AV*f{T) = Vffifa, T) - Vjjf(4> t , T) , 



rRen 



Rent 



(40) 



where 



V e f J n (<P,T) = V(<P)+T 



d 3 k 
J2nf 



In 1 - e 



-f3yJP+m?{<t>) 



(41) 



is the renormalized effective potential at finite temperature. 

Expression (38) for <r(T) can be easily identified as being the surface tension at finite 
temperature for the bubble wall. If one writes (38) as 



a(T) 



AttR 2 



[T eff ((j) wa i h T) -T eff ((j) f ,T)] 



(42) 
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where r e //(0, T) is given by 



r@ r ( - 1 + °° /■ d 3 fc - 1 

T e// (0, T) = ^ rfr J d\ \CeuM) + Yp E 7 (2^)3 ln K + + } ' < 43 ) 

T e ff((f>,T), as given by (43), is the effective action at finite temperature for a field config- 
uration at 1-loop order' 12 '. cr(T), given by (42), is the definition of the surface tension 
in field theory' 12 ^ . 

It is easily shown' 7 ! that in a high temperature limit ((3 — > 0, T — > oo) the terms in 
the expansion (33) that have higher superficial degree of divergence will contribute with a 
leading power in T. This contribution is just the first graph in (33) (m = 1 in (34)) such 
that at high temperatures one can consider only the following simple term of the graphic 
expansion in (33): 



r +°° r rf3 u i 

= / <?x [V"{4> b ) - V»m £ / ( " 2 ^r 2 | ^ > ( 44 ) 
J n=-oo J {^y' ul + k 2 + m 2 (<p f ) 

making the computation of AF(T) and then of T very simple with a bubble field confi- 
guration (f>b as given by (10). 

4 A Scalar Field cj) Coupled to Fermion Fields 

Let us now study a model of a scalar field coupled with massless fermion fields ip, with 
Lagrangian density given by 

C = i(^0) 2 - 1/(0) - + i$ffl>- (45) 

V(<p) in (45) is such that the finite temperature effective potential of the scalar field 
admits a first-order phase transition. For example with V((f>) given by Eq. (67), Sec. V. 
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The partition functions Z(4>b) and Z(cf>f) in (6) are now defined by the field expansions 
4>(x,t) — > <f>b(x) +r)(x,t) and ip(x,t) — > ip(x,t) for Z(0&) and (f>(x,t) — > 0/ + C(x, i) an d 
ip(x,t) — > ip(x,t) for Z((j)f). Again, 0j( t ) = 0/(t)(T) are determined by the minima of the 
effective potential V e ff((f>,T). Up to 1-loop order one gets the following expressions for 
Z(0 b ) and Z((j)f), respectively, 

Z(0 b ) = e- s ^ i ^^Dr ] exp|-^rfr|d 3 x^[-D £uc/ + y ,, (0 fc )]r / | x 

y DtfjDifj exp l-J^dr J d 3 xip [-ft- ig<p h ] ip J (46) 



x 



and 



Z(^) = e-^^^DCexpj-^drl^Chn^ + ^'^/Klx 



x j> DifjDipexp 



dr d xip[— ft — ig<j>f\ i/j 



(47) 



where the functional integration over the bosonic (fermionic) fields satisfy periodic (an- 
tiperiodic) boundary conditions in Euclidean time. 

As in (14) the free energy of the system can be written now as 



T — —T 



det(-n Eud + V"^ f )), 



det(- ft- ig<pb)(3 



det(- 



-AS 



(48) 



det(-D Eud + V"((/) b ))^ 
where AS = S Eucl (4> b ) - S Eucl {4> f ), with, from (45), S Eucl (<f>) = fg J d 3 x [\{d^f + V(</>) 

The determinant ratio of the bosonic part in (48) can be computed in exactly same 
way as before and we get the following expression for the nucleation rate T for the model 
(45), 



T = AT 4 



det(— ft — igcj) 



b)0 



_det(- ft-ig(j) f )fi_ 



exp 



AF(T) 
f 



(49) 
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where A and AF(T) have analogous identification with the terms appearing in (25). 
AF(T) can be given by both expressions (27) or (36) obtained by the eigenvalue method 
or by the graphic expansion, respectively. 

The fermionic determinant ratio appearing in (49) can be calculated by the same way 
as the bosonic one. We shall see that it provides also a correction to the exponential term 
of the nucleation rate T, that is, a fermionic correction to the factor AF(T) in eq. (49). 

If one uses the identity (which follows from charge-conjugation invariance): 



[det(- $ - igcf))] 2 = det(- # - ig<p).det(- $ + ig<j>) 



det 



(50) 



where l^xi is the 4x4 unit matrix. From (50) the denominator of the fermionic deter- 
minant ratio in (49) can be written as (0/ is the metastable vacuum field configuration) 



det(- - ig<j> f )p = [det(-a Eucl + g 2 4>))p\ 1 ■ (51) 

For the determinant involving the bubble field configuration one can again use a configu- 
ration like (10) for </> b , which describes a radially symmetric bubble solution. Making the 
following choice for the Dirac matrix in the radial direction 



7r = I 



12x2 
— 12x2 

where 12x2 denotes a 2 x 2 unit matrix, one can write det(— $ — ig^b(f))^ as 



V 



(52) 



where 



det(- p-ig<k( r )h = det fi (+) (0 fe ). det fi^^^) , 
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(53) 



& ± \ ( p b ) = -n Eud +g 2 ( pl±g 



d<j)b 
dr 



(54) 



The fermionic determinant ratio in (49) can be written then as (using again that 
IndetM = Tr In M ) 



det(- $ - ig^bh 
det(- $-ig<f)f)p 



= exp {Tr ln(- d - igfo)^ - Tr ln(- - ig<f> f )p} = 



= exp | Tr In 



- 2Tr In 



- a Eud + g 2 4>l + g- 



-n E uci + g 2 4>) 



dr 



+ Tr In 



n 2 ,2 



As in (31) one can write (55) in the form 



det(- - ig<j>b)p 
det(- $-ig<j>f)p 



exp < Tr In 



i + s p (4> f ) 



g\Vb - 0/) + g- 



dr 



+ Tr In 



1 + SMf) 



dr 



where 



+ 



(56) 



(57) 



- a Eud + m 2 F ((f) f ) 

is the fermionic analogous of Gp(4>f) given by (32) and m F (<j)f) = g 2 (p 2 . The expression 
in the exponent in (56) can be written as a series expansion like (34): 



Tr In 



Ati-<Pf)±g- 



dr 



+0O 



\m+l 



X 



E 

m=l 

+oo 

E 



Art - ti) ± g- 

1 



<9r 



x 



(58) 
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where u n = £^±^ ; for fermionic fields. As before, (58) can be expressed as a graphic 
expansion similar to (33) with the propagators Gp(4>f) exchanged now by Sp{(f>f) and the 
external lines given by g\<g - (j>)) + or g 2 ((f) 2 b - 4> 2 ) - g^ . 

When (56) is substituted in (49) and using the bubble field configuration (j) b as given 
by (10), in the thin- wall approximation one obtains 



T = ^T 4 exp 

where AF B+F (T) is given by 



AF B+F {T) 
f 



(59) 



AF B+F (T) = -^^AV B + F {T) + AitR 2 a B+F {T) , (60) 
3 

where AV e B f j F {T) = V B ^ F {(f) f ,T) - V B f j F {4> u T) is the effective potential difference be- 
tween the false and true vacua for the model (45). Up to 1-loop order one obtains, 



d 3 k 

4T J In (l + e^V^W)^ , ( 6 1) 



VSf^T) = V{^)+T f ^^{l-e^^^- 



d^k 

where m 2 B {(f)) = V"(<f>) and m 2 F ((f)) = g 2 (/) 2 are the boson and fermion effective masses, 
respectively, in the background field (f). a B+F (T) in (60) is given by 



- B+F (T) = [rf//(0^, t) - rf//(0 /5 d] (62) 



with 



rf//(0,T) = [drjd^- 



1 +oo 

(d,<p) 2 + v(<p) + - £ 

Z P n=-cx 



d 3 k 



In 



ul + k 2 + m%(<f>) 



- y 

y n=—cx 



d 3 k 



In 



ui + k' + m^+glVcPl 



(63) 
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T B ^ F (<fi, T) is the effective action, up to 1-loop order, for the model (45). As in (42), (62) 
defines the surface tension for the bubble wall. 

The determinat ratio in (55) can also be derived by the eigenvalue method if one knows 
how to compute the eigenvalues of the differential equation for the field configuration <fi wa ii 
describing the bubble wall, 



i-i 2/2 i dtPwall 

-o E ud + g <t> wa ii ± g dr 



(±) 0) 



Vwall 



,(±) 
'wall 



(j)]Viti(i) , 



(64) 



where, in momentum space, [e^//(j)] 2 = w\ + [E^ H (j)\ 2 . For the vacuum fields (fij and 
(fit we have the continuum eigenvalues (in momentum space) 



[e F {k)f = ul + [E F {k)f , [Ef (k)f = k 2 + 9 2 4>) 



(65) 



[e F {k)f = u? n + [E?{k)Y , [Ef(k)]* = k? +! W t 

for the false and true vacua, respectively. 

Using these eigenvalues in (55) one obtains an expression for (49) like (59), with 
AF B+F (T) given again by (60) and as it was shown in the bosonic case, from (65), the 
expression for AV e B J F (T) remains the same as that one given in (60). The expression for 
a B+F (T) is now rewritten in terms of the eigenvalues for the bubble wall field configuration 
<fiwaih such that, from (64) and taking into account the bosonic part of a B+F {T) given by 
(29), one can write 



r B+F 



CO = ^( T )-Z^i E ^ri(j( ± )) + 21n(l + e- 

f U(+) >•?'(-) 



]T[2/3£f(;)+4m (l + e -^/«) 



(66) 



where <r B (T) is given by (29) and E F (i) = \jk 2 + m 2 F (<fif). 
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5 Simple Example 



Let us now illustrate our method for a simple model that exhibits metastability and for 
which we can estimate the eigenvalues related with a bubble wall field configuration 4> wa ii- 
Consider the model of a scalar field 4> coupled to fermion fields ip with Lagrangian 
density given by (45) and potential V(<f>) given by 

v(<f>,h) = ~^ + ^-h < l>, (67) 

where h is assumed to be a constant external current. 

It is well known that the potential (67) exhibits metastability when one varies the 
sign of /it 3 ]'' 13 ]. V((f),h) has two minima, <f>f(h) and t (/i). For h > (j) t describes the 
stable phase and <fij the metastable phase. For h < the rules of <f>f(h) and (f>t(h) are 
reversed. At finite temperatures, we can always find values for the Yukawa coupling g and 
for the constant h, such that the condition of existence of metastability is satisfied' 14 ] , with 
(f>f(t)(h,T) describing the metastable (stable) phases, respectively. In the limit h — > ± , 
the vacua solutions (f>f,t(h,T) tend to the limit ±<p ( T ) = ± (^r^)", where m 2 (T) is 
the finite temperature effective mass for the scalar field <fi ( m the high-temperature limit, 
m 2 (T) ~ m 2 — — 2 |-. ) and one of the non-trivial static solutions of the Euler-Lagrange 
equation for the field in model (45), with potential given by (67), is 

(m(T) \ 

4>kink{x L ) = o (T)tanh I — j^^l 1 , (68) 

which describes a planar interface between the two vacua ±4> (T) (h — 0). xl is the 
longitudinal component of the spacial vector ( X — (^Li ^t) )• 

For suficiently small values of h, the bubble field configuration 0&(r) is well approxi- 
mated by the radially symmetric solution f 13 ': 
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1 1 
M r ) - + W + 2^f- W tanh 



m(T) 
"7T 



R) 



(69) 



which describes a bubble of radius R of the nucleating vacuum (f> t = 0t(T) (the true 
vacuum) embedded in the metastable vacuum </>/ = 4>f(T) (the false vacuum). From (69), 
the thin-wall approximation is equivalent to considering R » m _1 (T), where Ai?, the 
bubble wall thickness, is proportional to m _1 (T). 

From (29) and (66) we see that we have to compute the eigenvalues of the differential 
equations (approximating for small values of h) 



V 2 + g 2 (t? wall ±g\V<p wa u\ 



B __ rp2 B 
rwall wall t wall 



(70) 



Twall — 



E 



(±) 

wall 



Vwall 



for the bosonic and fermionic terms, respectively. 

If one writes Lp wa ii = ^ n ,i(r)&i, m (0, ip) , and as (p wa ii can be written as a radially 
symmetric solution, (70) is given by 



d 2 



2d 1(1 + 1) A 



dr 2 r dr 



Qm 



tf B| ,(r) = El^^ir) 



(71) 



d 2 2 d 1(1 + 1) 2 2 d(/w/ 

-— + n + 9 <Pwall ± # 



= 



E 



(±) 



*5 } (r) 



o?r 2 r r 2 wa dr 

for the bosonic and fermionic terms, respectively. 

For small enough h one can consider the bubble wall field configuration 4> wa u as given 
basically by a field configuration like (68), that is, a kink-like field configuration and (71) 
is reduced to 
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ar 2 r ar r 2 



dr 2 r 



2 d + m 2 (T) 



5 2 - 5(5 ±l)sech 2 



V2 



' m(T) 

~7T 



1 7 



<i*»,i(0 



(72) 



(±) 



n,Z 



*5 } (r). 



In (72) 5 2 = As 5 and A are arbitrary constants in the model, 5 is an arbitrary 

positive constant. Choosing 5 as a positive integer number, the eigenvalues of (72), for low 
temperatures, can be easily estimated (the differential equations in (72) are equivalent to 



the Schrodinger equation for a Posch- Teller potential). The eigenvalues E 2 x and 
are given by (see ref. [15] for instance) 



E 



(±) 



E l,i 



l(l+l)-2 
R 2 ' 



n 



n 




1 



(73) 



k 2 + 2m 2 + + n -> fc (continuum) 



and 



(±) 



(74) 



./.j.^ , ^ 2 0,1,2, ...,5-1, for (+) 

^ + mM 5-f) + 0(^), P = 

[ 1,2, ...,5-1, for (-) 

k 2 + ^ + §m 2 5 2 + O(^) 

where, in the equations above, i? is the bubble radius. 

For the bosonic eigenvalues we have e 2 os = (^p) + -E" 2 ^ (at T ^ 0). Therefore , 
for j — 0, £"q corresponds to a negative mode, which is associated with the instability 
of the critical bubble. For j = 0, Eft-y (m = 0, ±1) represents the three translational 
modes related to the space translation invariance of the bubble. At high temperatures 
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we still may associate the negative and the zero eigenvalues with l ^ l+ ^~ 2 , for I — and 
I = 1, respectively. Note that, at T = T c , the critical radius of the bubble must go to 
infinity, R cr (T c ) — > oo, and therefore the negative eigenvalue E% = E 2 _ = — vanishes, 
as expected. 

Taking the limit — > ± is equivalent to taking the limit of infinite bubble radius, 
R — > oo. In this case the bubble wall solution <fi wa u tends to the planar wall interface 
solution (ftkink, eq. (68), and the eigenvalues (73) and (74) become that of the planar wall 
solution. (In (73) and (74), in the limit /z — ^ ± (R — > oo), the relevant contribution of 
the /-dependence comes from large values of / such that the sum over /, of the eigenvalues, 
J2i^o (2Z + 1) , can be replaced by a continuum integration in the two-dimensional momenta 
parallel to the wall's surface, / cPk^. ) 

The evaluation of expression (62), for the surface tension, in terms of the eigenvalues, 
of course, is not an easily task. However, as shown in the last two sections, we still can 
use the expansion for the determinants, Eqs. (34) and (58). For example, in the limit 
h — > 0^, we can easily estimate the extra contribution to the surface tension, coming 
from fluctuations around the wall field configuration (f) wa ii- From (34) and (58) one gets, 
in leading order in T (that is, from the tadpole graphs, like the one in (44) ), that these 
extra contribution for the surface tension cr , given by (18), is given by (approximating 
Kaii{r) by (68), h -> 0±) 

0"r = &t + a T 5 (75) 

where 
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l r+°° . r „ r c? 3 fc 



dx L [V"((f) kink (x L )) - V"(4>q)} 



(76) 



for the bosonic contribution and 



/+oo 
d,XL 
-co 



d 3 k 



^ I (2^ 



\Jk 2 + m 2 F ((/) ) (e^V^+^F^o) 



+ 1 



is the fermionic contribution. In (76) and (77), <pkink(xL) is given by (68), 4>q 

mU<M = ^U=MT) and m^o) = 9 2 4>l(T) = ^f^" 

In the high-temperature limit, T >> m (A ~ £/ 2 , A << 1), we obtain that 



(77) 



6m. 2 (T) 
A ' 



a T ~ 



V2m(T)T 2 



+ 0(T) , 



<r£ ~ + C(T) . 



A 



Using that S" 2 = we obtain 



(78) 



0T — 



V2m(T)T 2 



q2 

1 + T 



+ C(T) 



where m 2 (T), in leading order in T, is given by m 2 (T) ~ m 2 — ^ — 



(79) 



6 Conclusions 

In this paper we have shown how to compute the nucleation rate in field theory at finite 
temperature and also we have shown how the prefactor term (the ratio of the determinant) 
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contribute to the exponential term in the classical nucleation rate, showing that it provides 
a finite temperature quantum correction to the exponential term. 

The procedure we have used is completely consistent and it is based on the usual 
description of the decay of metastable states developed first by L anger®, translated here 
for finite temperature field theories. Actually, the final expressions that we have obtained 
for the nucleation rate, Eq. (25) or Eq. (59), are similar the one proposed by Langer, 
in which the droplet energy AE is replaced by the droplet free energy AF(T), when one 
takes into account the prefactor term. 

For the evaluation of the prefactor term (the determinant ratio) we have used two 
methods. The first method is based on the explicit evaluation of the eigenvalues. The 
second procedure is based on a field theoretic expansion for the determinant ratio, such 
that we avoid the task of computing the eigenvalues of the differential equations. Both 
methods depend on the knowledge, at least in an approximate way, of a field configuration 
describing the bubble wall. 

We have shown (using the thin-wall approximation) that the nucleation rate has a 
form that is much like the "classical" one (without taking into account the determinants). 
In the general expression for the nucleation rate T, we have shown that AF(T) is given 
by 

AF(T) = -—R 3 AV eff (T) +AnR 2 a(T) , 

where AV e ff(T) is just the effective potential contribution that contains the "classical" 
result plus a quantum correction that is calculable in field theory. The surface term cr(T) 
is a sum of two contributions: a "classical" term plus a temperature dependent quantum 
correction. The general structure of o(T) is given by Eq. (66). 

Due to the similarity between AF(T) and the free energy of the bubble, there is appar- 
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ently no distintion between our results and analogous ones presented in the literature' 1 ''' 5 '. 
This will be true as long as one computes AV e ff(T) and cr(T) in a consistent way. We 
have given in this paper examples of how to carry on such a calculation consistently. 

It is interesting to contrast the expression that we get for the nucleation rate with 
expression (1). We have shown that all quantum corrections, due to fluctuations around 
the bubble solution, to AS* in the exponential term in (14) are in the determinant ratio. 
For the nucleation rate associated to a critical bubble (when AF(T) assumes its maximum 
value), from (25) and (26) we get the following complete expression for T: 



(T) 



Io(R cr )e ""'^ 



(80) 



with the preexponential term I (R cr ) given by 



lo(Rcr) 



1 



V2TR cr (T) 



TTsin 



AE(R a 
2ttT 



.V2TR cr (T) 

In (80) we have used (from the expression for AF(T)) that 



rp4 



(81) 



R cr (T) 



2a(T) 



and in (81) we have used, from Eq. (73), that the negative eigenvalue is given by 



E 2 



Rcr(T) 



2, for the model (45) with potential given by (67). 



We remark also that in ref. [6] the preexponential factor Io(R) for critical bubbles was 
estimated, in the case that T represents the probability per unit volume per unit time to 
nucleate a critical bubble of the hadronic phase inside the quark phase, to be' 6 ' 



where ( q and rj q are viscosity coefficients, £ q is the correlation length in the quark-gluon 
phase, Auj is the difference of enthalpy density between the quark-gluon phase and a is 
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the surface tension, assumed there to be temperature independent. The authors in [6] 
have pointed out then that the usual approximation for the preexponential term, like 
T 4 or T c 4 , can not be a good approximation, since, from (82), the preexponential is a 
complicate function of the temperature, diverging as T — > T c . The same is true here, 
where the preexponential factor is given by Eq. (81). 

Despite the fact that we have carried out our analysis by using a thin-wall approxima- 
tion, that is only valid for nucleation happening at temperatures close to T c , or equiva- 
lently, for small supercooling, it is possible to extent the method for other cases of interest, 
for example, thick-wall bubbles' 5 ' 16 ', provided that we can find an approximate solution 
4>b(f), describing the bubble field configuration. 
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